Quantum Stress: Density Functional Theory Formulation and Physical Manifestation 
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The concept of "quantum stress (QS)" is introduced and formulated within density functional 
theory (DFT), to elucidate extrinsic electronic effects on the stress state of solids and thin films in 
the absence of lattice strain. A formal expression of QS (cr*^) is derived in relation to deformation 
potential of electronic states (H) and variation of electron density (An), cr'^ — HAn, as a quantum 
analog of classical Hook's law. Two distinct QS manifestations are demonstrated quantitatively by 
DFT calculations: (1) in the form of bulk stress induced by charge carriers; and (2) in the form of 
surface stress induced by quantum confinement. Implications of QS in some physical phenomena 
are discussed to underlie its importance. 

PACS numbers: 62.20.-x, 75.15.-m, 64.60.-i, 68.35.Md 



A fundamental property of solid materials is their 
stress state. At the equilibrium lattice constant, the bulk 
of a crystalline solid is stress free, but the surface of a 
solid has intrinsic non-zero stress and stress is commonly 
induced by any form of lattice distortion [l[ . The stress 
(strain) state of a solid or thin-film material has profound 
effects on its thermodynamic stability and physical and 
chemical properties [l|-|3|, and has been employed in a 
wide range of applications such as electromechnical de- 
vices Q, mechanochemcial sensors Q and flexible elec- 
tronics p , and even to make new nanotructures 0, @| . 
Here, we introduce the concept of " quantum stress (QS)" , 
which adds an interesting quantum mechanical aspect to 
our conventional view of "classical stress (CS)" based on 
classical mechanics. We mathematically formulate the 
expression of QS within density function theory (DFT), 
and use DFT calculations to demonstrate quantitatively 
two distinct physical manifestations of QS, in the form of 
bulk stress induced by charge carriers in a homogeneous 
system of crystalline solids and in the form of surface 
stress induced by quantum confinement in a heteroge- 
neous system of nanoscale thin films. We will then apply 
the concept of QS to elucidate a few examples of physical 
phenomena that underlie the importance and usefulness 
of QS. 

Concept of QS. We first introduce the concept of QS 
in contrast with the CS. Figure 1 illustrates the funda- 
mental difference between the QS and CS using a sim- 
ple model of a one-dimensional (ID) lattice. Consider 
a lattice is under compressive (Fig. la) or tensile lat- 
tice strain (e), such as in an epitaxial film due to lattice 
mismatch between the film and substrate U,^. The 
" atomic" lattice deformation energy can be expressed as 
E — {l/2)Ye'^V, where Y is Young's modulus and V is 
the volume of lattice. By definition, the lattice forma- 
tion induced lattice stress, which we refer to here as CS, 
is expressed as a'^ = {l/V){dE/de) — Ye , the Hook's 
law. Now, consider an equilibrium lattice in the absence 
of strain (e = 0), but electronically perturbed or excited, 
such as an electron is kicked out by a photon leaving be- 



hind a hole, as shown in Fig. lb, which redistributes the 
electron density. The "electronic deformation" energy 
can be expressed as E — fJ-AN, where /i is electron chem- 
ical potential and AN is the change of number of elec- 
trons. Then, the lattice stress induced by the electronic 
deformation, which we refer to as QS, can be expressed 
as = {l/V){dE/de) = EAn, where S = dji/de is de- 
formation potential and An is change of electron density. 
The expression of cr'^ = SAn can be viewed as a quan- 
tum analog of Hook's law. Below, we provide a formal 
derivation of CS and QS within DFT. 

DFT formulation. Following DFT [lH, the total en- 
ergy functional of a solid is written as 

S[n(f),{i?,„}] = E,[n{r)]+E,,t[n{r)ARn^}] + Ei[{R,n}] 

(1) 

Ee[n{r)]is the electronic energy functional of charge 
density n(r, including kinetic and electron-electron inter- 
action energy, Eext[n{r), {Rm}] is the ion-electron inter- 
action energy, -E/[{-Rm}] is the ion-ion interaction energy 
and {Rm\ are atomic coordinates. First, for complete- 
ness, let's rederive the Hook's law, the expression of CS 
under finite lattice strain. Consider coordinate transfor- 
mation r = (1 -f {eij})'r^ and Rm = (1 + {£ij})Rm under 
strain {Sy } , where and {-R^} are the electronic and 
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FIG. 1: (Color online) Schematic illustration of CS versus QS. 
(a) The CS (f"^) induced by applying a compressive lattice 
strain (e). Arrows indicate stress and force directions, (b) 
The QS (u'^) induced by a hole excited by a photon. 
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atomic coordinates of strain-free equilibrium lattice. Let 
n°(r^) and be the ground-state charge density be- 

fore and after strain is applied. By definition, the stress 
tensor is expressed as 



C _ 1 dE[n{r),{R„,}] 
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where Er = E^xt + Ej. Since n^{r) is the ground-state 
electron density at r and {Rm} , according to Hohenberg- 
Kohn theorem [ll|, we have 
and Eq. ([2]) becomes 
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For simplicity, assuming hydrostatic strain, £ij — eSij , we 
expand Er in e 

ER[n^, {R,n}] - ER[n^, {R^J] + e V i?," ' " " 
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Then the stress can be expressed in the first order of e 



as 



where Y 
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(5) 

is the Young's 



modulus. 

Next, we derive the QS induced by electronic exci- 
tation and perturbation without applying lattice strain 
{sij = 0). Consider a variation of electron density 
from rfi{f^) the ground-state density at f* and {R^} 
as n*{f^) — n"(7*) -|- 6n{f^) . (Below, for convenience, 
we will neglect the superscript for r*.) The differentials 
of energy functionals are 

FK(f)]=F[nO(f)]+ / f^^i^) &n{T)df, (6) 
Jv V 5n{r) y„o 

The stress tensor is 

Q _ 1 dE[n{r),{Rra}]^ 
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Where = 9„(i?e + Ef-xt) is the electron chemical po- 
tential. To arrive at Eq. ([7]), we used the condition 



that the strain- free ground-state solid is stress free, i.e., 
( 1~ ) = . It can be shown that the second term 

in Eq. ^ vanishes because chemical potential remains 
uniform and the number of electrons is independent of 
strain, so we have the expression of QS as 
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In a homogeneous crystalline solid, to a good approxima- 
tion, the electron deformation potential S = d^/dsij is 
uniform as the electron density remains uniform before 
and after strain is applied. Then, the expression of QS 
can be simplified as 



EAn 



(9) 



Equation ^ can be viewed as a quantum analog of Eq. 
([5]), with , S and An playing the role of cr'^, Y and e, 
respectively. However, Eq. ([8|) must be used if S is not 
uniform in a heterogeneous system. For example, in thin 
films (hetero-junctions) when strain is applied, charge 
will be redistributed in the surface (interface) regions due 
to the nonuniform 

Notice that the CS reflects the atomic and lattice size 
effects. Although in the derivation of CS (Eq. ([5])), 
the energy expression of Er has implicitly the electronic 
contribution through the ion-electron interaction energy 
Eext as strain changes electron density from n°(f*) to 
n(r) , the stress, nevertheless, follows the Hook's law 
of classical mechanics depending only on atomic coordi- 
nates. In other words, even though the atomic and lattice 
size is associated with electronic structure, the effects of 
the ground-state electronic structure can be cast into the 
atomic and lattice size effect, having a CS manifestation. 
It is for this reason that the CS can be modeled by em- 
pirical interatomic potential involving explicitly only the 
atomic degrees of freedom as done in molecular dynam- 
ics simulations. In contrast, the QS has a pure electronic 
origin involving explicitly only the electronic degrees of 
freedom [n(r)] that cannot be cast into the atomic or lat- 
tice size effect. Consequently, the QS must be described 
by the quantum mechanics of the electronic structure 
rather than the classical mechanics of the atomic struc- 
ture. Below, usin g fi rst-principles quantum-mechanical 
stress calculations [13, we directly quantify the magni- 
tude and reveal the nature of QS in two distinct mani- 
festations. 

QS induced by charge carrier. We first demonstrate 
the QS for the case of a homogeneous system where 
Eq. ^ can be applied, in the form of bulk stress when 
an electron is added to or removed from a solid lattice, 
as in the case of semiconductor doping or photo-excited 
charge carriers in solids. We have performed DFT calcu- 
lations il3:| of the lattice stress induced by adding elec- 
trons and/or holes (i.e., removing an electron) to a fi- 
nite lattice of Al (metal). Si (elemental semiconductor). 
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FIG. 2: (Color online) The QS induced by electrons (trian- 
gles) and lioles(circles) as a function of carrier density in (a) 
Al, (b) Si, (c) GaAs, and (d) Zr02. Solid lines are linear fits 
to the data, and the dashed lines are extensions of the solid 
lines to indicate different slopes for electrons from holes. 



GaAs (compound semiconductor) and Z1O2 (insulator), 
and graphite with a hexagonal lattice. Figure 2 shows 
the calculated ct'^ as a function of An for Al, Si, GaAs 
and Zr02, which shows an almost perfect linear depen- 
dence for all the cases, in excellent agreement with our 
theoretical derivation of Eq. ^ . In general, electrons in- 
duce compressive QS (by convention, compressive stress 
is defined as negative); while holes induces tensile QS. In 
plotting Fig. 2, the QS values are taken from the diag- 
onal terms of the stress tensor along the principal axes, 
since stress is isotropic in a cubic lattice. More generally, 
electrons or holes may induce anisotropic stress, such as 
in a hexagonal lattice of graphite (see Fig. 4 below). 

According to Eq. ([9]), the slope of cr'^ vs. An equals 
to the deformation potential, S. For a metal, S = 
is the same for electron and hole because of the electron- 
hole (e-h) symmetry in the metal, as seen for Al in Fig. 
2a, and we found Eai = —10A9eV. For a semiconduc- 
tor or insulator, however, the deformation potential for 
electron (S*^ = ^^^^m ^ CBM stands for conduction band 



de 



minimum) is different from that for hole (S 
VBM stands for valence band maximum) because of 
the e-h asymmetry, as seen for Si, GaAs and Zr02 m 
Fig. 2. We obtained that = -8.65, = -9.51; 
^GaAs = —9.77, =^GaAs = —7.33; ^%r02 ~ —12.36, 
^Zr02 ~ ~8.87, which are in good agreement with pre- 
vious results [IJ]. In general, the larger the band gap is, 
the larger the e-h asymmetry and hence the larger the 
difference between the electron and hole deformation po- 
tential will be. For Si, S*^ is smaller than S'', possibly 
because it is an indirect semiconductor, while for GaAs 
and Zr02, S"^ is larger than S''. 

We note that conventionally, the deformation poten- 
tial is derived by calculating the valence and conduction 
band edge positions as a function of strain, which can be 



difficult for DFT methods because of the arbitrariness in 
the absolute value of band energy. Here, our QS calcula- 
tion provides an efficient and effective method to derive 
the deformation potential from total energy calculations 
without the need of calculating band structure. 

QS induced by quantum confinement. We next demon- 
strate the QS for the case of a heterogeneous system 
where Eq. ([8|) must be applied, in the form of surface 
stress of nanostructures, as in the case of epitaxial growth 
of nanometer-thick metal thin films when the quantum 
size effect is prominent. All the crystalline solid surfaces 
have a non-zero intrinsic surface stress [ll,[3]. Classically, 
one may view this by imagining that the surface layer 
of atoms would have a smaller lattice constant (bond 
length) if they are separated by themselves; they have 
to be stretched apart when placed on top of the film to 
match the lattice constant of underlying film. For this 
reason, most solid surfaces have a tensile surface stress 
with a well-defined magnitude, a characteristic surface 
property of a given film structure. However, if the thick- 
ness of a film is reduced to nanoscale comparable to the 
de Broglie wavelength of electrons, quantum confinement 
becomes prominent giving rise to the formation of dis- 
crete quantum well states, known as the quantum size 
effect (QSE) [H, [i3l_ The QSE has been shown to mod- 
ify surface energv~[l8| . Here, we demonstrate that QSE 
will also modify surface stress, as a distinct manifestation 
of QS induced by quantum confinement [l9l |. 

Figure 3 shows the calculated surface energy (7) and 
surface stress (a) as a function of Pb(lll) film thick- 
ness (d). 7 displays an oscillatory dependence on d, as 
known before [18|. What's new is that a displays a much 
stronger oscillatory dependence on d, showing a much 
stronger QSE on surface stress than on surface energy. 
The thickness dependence of the quantum surface stress 
can be understood from the thickness dependence of the 
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FIG. 3: (Color online) The calculated surface energy and sur- 
face stress of Pb(lll) film as a function of film thickness, 
demonstrating the surface QS in metal nanofilms induced by 
quantum confinement. 
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quantum well states formed in the thin film, which mod- 
ulates the thin film deformation potential (S) and surface 
charge density (An) as a function of the film thickness, 
leading to the oscillating quantum stress. Because both 
S and An are nonuniform in thin film, the simplified ex- 
pression of Eq. © cannot be used. The results in Fig. 3 
are the integrated results of Eq. (|8]) for each film thick- 
ness. Empirically, however, we may conveniently divide 
the surface energy into classical and quantum contribu- 
tions as 7 = -I- (d) . Then, by definition, we express 
surface stress as 



f ^7 _ 1 d-y'^ I d-y'^ 
Ade A de A de 



■id) (10) 



where ^is surface area, which is also divided into clas- 
sical (fT*") and quantum contribution (cr'^). 7"-^ and a*^ 
represent, respectively, the classical surface energy (bond 
breaking energy) and surface stress (bond deformation 
energy) of a macroscopic thick film independent of film 
thickness; 7*5 and cr'^ represent, respectively, the quan- 
tum surface energy and surface stress, arising from quan- 
tum confinement in a nanoscale thin film, as a function 
of film thickness d. As the film thickness increases, 7*^ 
and cr*^ will eventually diminish and the system resumes 
the classical behavior. 

Implications of QS. The concept of QS has implications 
in a broad range of physical phenomena and technologi- 
cal applications that are based on coupling of electronic 
structure with lattice strain. We have already shown 
that the DFT calculation of QS provides an effective 
method for deriving deformation potential without the 
need of calculating band structure, which circumvents 
the difficulties encountered by previous methods as well 
as saves computational time. Physically, the QS induced 
by charge carriers will help us to better understand the 
physical nature of semiconductor doping in terms of the 
dopant-induced lattice stress, by differentiating the QS 
induced by electrons and holes from the CS induced by 
size difference between dopant and host atoms In 
general, it is easier to dope an element whose QS and 
CS compensate each other, i.e, small n-type dopants or 
large p-type dopants, which induce the smallest overall 
amount of stress. There is a partial cancellation effect on 
the QS between the n- and p-type dopant, which makes 
the co-doping of both types of dopants easier. 

Another QS-related physical phenomenon is the pulse 
laser induced structural phase transition, such as the 
graphite-to-diamond transition 2l!, i223. In a pulse lasing 
experiment, a high density of charge carriers (electrons, 
holes and excitons) is photo-excited in a small volume 
for a very short time. We argue that such charge carriers 
exert a large QS to the local lattice, causing effectively a 
"pressure-induced" structural phase transition. To sup- 
port our point of view, we calculated the QS exerted by 
the photo-excited carriers (holes in the valence bands plus 
"free" electrons) to an ABC-stack or AB-stack graphite 



lattice [12[. The QS is found to be tensile and highly 
anisotropic with the largest component along the z-axis 
and all the components increase approximately linearly 
with carrier density (Fig. 4a for the ABS-stack graphite). 
This is because the QS is dominated by the contribution 
from holes in the valence band of Pz orbital, as indi- 
cated in Fig. 4b. The magnitude of the QS induced by 
a single hole in the 6-atom cell is as high as 20-30 GPa 
(Fig. 4a), which indicates that the pulse laser can induce 
a huge "local" stress (pressure) in the graphite lattice, 
larger than the critical pressure needed for the graphite- 
to-diamond transitions [23j. Furthermore, we relaxed the 
graphite structure under the QS exerted by the charge 
carrier, we directly observed the lattice transformation 
of graphite into cubic diamond as the QS is gradually 
decreased upon structural optimization [12| . These re- 
sults shed new lights on the understanding of the pulse 
laser induced graphite-to-diamond transition. More gen- 
erally, we expect that the QS generated by the excited 
carrier is a key physical parameter in understanding a 
wide range of phase transitions induced by radiation of 



energetic particles [24 



The realization of the surface QS in nanoscale thin 
films has a profound effect on film stability and growth, 
especially when the film is strained. The QSE has been 
shown to be dominant in the gro wth of metal thin films 
on semiconductor substrates |25|. On the other hand, 
the strain effect is ubiquitous in heteroepitaxial growth 
So far, however, these two important effects have 
been usually studied separately focusing on one while ne- 
glecting the other, despite the fact they are both present 
in many cases. One outstanding difficulty is the lack of a 
theoretical framework underlying the fundamental rela- 
tionship between the two effects. The introduction of the 
surface QS establishes a direct link between the quantum 
size and strain effects on the surface energy of thin films, 
underlying the interplay between the quantum size and 
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FIG. 4: (Color online) (a) Stress components (the z- 
component along the c-axis) induced by charge carriers as 
a function of carrier density in an ABC-stack graphite lat- 
tice, (b) The hole partial density of states with one electron 
removed from a 6-atom hexagonal graphite unit cell. The in- 
sets show the corresponding differential hole charge density 
of Pz orbital shape. 
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strain effects on thin film stability and allowing a theoret- 
ical prediction of the thickness-dependent surface energy 
in a strained quantum film. 

In conclusion, we present a rigorous derivation of QS 
within DFT, which is shown in linear relation to deforma- 
tion potential and charge density variation, as a quantum 
analog of Hook's law of CS. Two distinct physical man- 
ifestations of QS induced by charge carriers and quan- 
tum confinement are demonstrated by direct DFT cal- 
culations, confirming the theoretical formulation. The 
concept of QS adds an important quantum mechanical 
aspect to our classical view of stress. It manifests itself 
broadly in various mechanical properties of different ma- 
terial systems, especially in low-dimensional nanostruc- 
tures where quantum effects become prominent [l^. It 
will be especially relevant in advancing our fundamen- 
tal understanding of physical phenomena and technolog- 
ical applications that are based on coupling of electronic 
structure with lattice stress, such as charge carriers con- 
fined in quantum wells, wires and dots, free carriers cre- 
ated in nanofilms by electrical gating, electroelastic ef- 
fects, magnetoelastic effects and biological cell deforma- 
tion due to charging and polarization. 
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